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Replica Limit of the Toda Lattice Equation 
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In a recent breakthrough Kanzieper showed that it is possible to obtain exact nonperturbative 
Random Matrix results from the replica limit of the corresponding Painleve equation. In this article 
we analyze the replica limit of the Toda lattice equation and obtain exact expressions for the two 
point function of the Gaussian Unitary Ensemble and the resolvent of the chiral Unitary Ensemble. 
In the latter case both the fully quenched and the partially quenched limit are considered. This 
derivation explains in a natural way the appearance of both compact and noncompact integrals, 
the hallmark of the supersymmetric method, in the replica limit of the expression for the resolvent. 
We also show that the supersymmetric partition function and the partition function with fermionic 
replicas are related through the Toda lattice equation. 



Introduction. The replica trick originally introduced 
in JjJ has been widely used in the theory of disordered 
systems ranging from spin glasses to QCD Q). The 

quest for its rigorous understanding has been ongoing 
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in this letter. 
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lj,|15[ and will be further pursued 
The idea is to simplify the calculation 
of the disorder average of the logarithm of the partition 
function starting from 



' Z n — 1 
<logZ)= lim ' ' 

n— »o \ n 



(1) 



This method should work if we know the right hand side 
of this equation as an analytical function of n. In prac- 
tice, though, one is usually only able to calculate (Z n ) 
for positive integer values of n which is not sufficient to 
reproduce the complete function. For example, terms of 
the form (sin7rnx)/n or singularities along the real n- 
axis could lead to a breakdown of the replica limit. The 
first problems with the replica trick appeared in the the- 
ory of spin glasses Q where the result for the entropy 
turned out to be not positive definite. These problems 
were resolved by means of an intricate scheme of replica 
symmetry breaking ||. Although the replica trick has 
been used widely in the theory of disordered systems, 
it has been criticized as well ||[|. I n particular, it has 
been a long standing belief that the replica trick can only 
be applied to asymptotic expansions which can be re- 
summed in some cases jn|. This belief was shattered 
in recent breakthrough by Kanzieper [l4| . In this letter 
we will obtain further exact nonperturbative results from 
the replica limit of the Toda lattice equation. Both this 
method and the method of |Q rely on the relations be- 
tween Random Matrix Theories and the theory of exactly 
integrable systems. 

We will analyze the replica limit of correlation func- 
tions of a Random Matrix partition function. In the case 
of the resolvent we start from the identity 



Tr- 



1 



D 



1 ,. Q /det n (D + x) 
= — lim d x ' 

N n->0 



Here, D is an N x N random matrix and the average is 
over the probability distribution of its matrix elements. 
The replica limit can be calculated in two different ways, 
by calculating the partition function for fermionic repli- 
cas, i.e. for positive n, or by using bosonic replicas, i.e. 
for negative n. 

An exact method that avoids difficulties of the replica 
trick is the supersymmetric method |jl6 17,ljjfl. In that 
case the resolvent follows from the identity 



G(x) 
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— ( 
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I det(Z) + x + J) 
' 7 \ det(D + x) 



(3) 



j=o 
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This method is not applicable, though, to cases where 
the partition function cannot be expressed as the average 
of a determinant, making it essential to understand the 
limitations of the replica limit. In this letter we will not 
discuss other alternatives to the replica method such as 
for example the dynamical method |l9fl . 

Both the generating functions for the replica trick and 
the supersymmetric method can be rewritten as a non- 
linear (T-model. Generally, it is believed HJ|^^,0,||| 
that the replica limit can only be applied to asymptotic 
expansions such as the expansion in 1/x of the resolvent. 
In H it was argued that the replica trick even fails to 
reproduce the oscillatory terms in the asymptotic expan- 
sions of correlation functions. Indeed, this is the case if 
one only takes into account the leading saddle point man- 
ifold. The oscillatory terms in the asymptotic expansion 
can be obtained from a subleading saddle point manifold 
H (see also @||). However, exact results could only be 
derived when the asymptotic expansion terminates . 
In H this was related to the Duistermaat-Heckman lo- 
calization property of the saddle point integrals. 

A real breakthrough was recently made by Kanzieper 
He succeeded to obtain exact nonperturbative re- 
sults by means of the replica method. His approach was 
based on the property that the generating function for 
the resolvent or the two-point correlation function can 
be expressed in terms of a solution of a Painleve equa- 
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tion where the replica index only occurs as a parameter 
in its coefficients. The replica limit of the resolvent is 
then obtained from solution of this Painleve equation for 
n — * Q. In particular, this allowed him to reproduce 
exact nonperturbative results with an asymptotic series 
that does not terminate. A first hint that it is possible 
to obtain such exact results surfaced in |Tl]] where the 
leading logarithmic singularity of the small x expansion 
of the resolvent was obtained from the replica trick. 

In this letter we will analyze the replica trick for the 
two-point function of the Gaussian Unitary Ensemble 
(GUE) and for the microscopic limit of the resolvent of 
the chiral Unitary Ensemble (chUE) — the random ma- 
trix ensemble relevant for QCD 

The two-point function of the GUE can be expressed 
in terms of a generating function for n fermionic replicas 



G(r) = - lim 



1 d 2 Z n (ir) 



n 2 dr 2 



(4) 



where the eigenvalue representation of the nonlinear a- 
model for the generating function is given by 0] 



/l n 
; J d\ k e 
- 1 k=\ 



-^A 2 ({\ k }). (5) 



Here, A({X k }) = Y\ k< i(X k — A;) is the Vandermonde de- 
terminant. For n bosonic replicas, the two-point func- 
tion is still given by (^) but with a generating function 
Z_ n (r), that now can be expressed as an integral over a 
noncompact manifold ||, 

ZMir) = Z^Yjf I UdX k e^A 2 ({X k }). (6) 



The sign of the exponent is determined by the sign of the 
infinitesimal increment of ir. The fermionic generating 
function satisfies the Toda lattice equation M 



Z 2 n (ir)d 2 log Z n (ir) 



(n + 1) 2 



(ir)Z n _ x (ir). (7) 



If we interpret the generating function for negative n as 
the the bosonic generating function (||) this equation is 
valid for all integer values of n. In the replica limit, the 
l.h.s. of this equation is n 2 times the connected two-point 
function. On the r.h.s. we find from (pf) and (H) that 



Zx{ir)Z^i{ir) = 2 



. sin r e v 



(8) 



which is the analytical result for the two-point function 
pj ! In H this result was obtained from the replica limit 
of the Painleve V equation. The present method nicely 
explains the factorization of the two-point function into 
a compact and a non-compact contribution. In Q it was 
argued that the replica limit reproduces the correct result 
in this case is because the asymptotic expansion in 1/r 



terminates. We will now analyze the microscopic limit of 
the resolvent G(x) of the chUE for which the asymptotic 
expansion in 1/x does not terminate. Nevertheless, we 
will obtain the exact analytical result for G{x) from the 
replica limit of a Toda lattice equation. 

For topological quantum number v the microscopic 
limit of the resolvent for the chUE is given by 

G(x) = x{K l/ {x)I u (x)+K v _ 1 (x)I u+1 (x)) + -. (9) 

x 

where K v (z) and I u {z) are modified Bessel functions. 
The most direct way to obtain this result is the super- 
symmetric method ^3[|4 15 where the appearance of 
the compact and noncompact integrals can be related to 
the symmetries of the underlying partition function [ pT| . 
However, it was first obtained in | po| by integrating the 
microscopic spectral density which can be easily derived 
from the orthogonal polynomial method ^,^,^8|. In 
this result was obtained by solving the replica limit 
of the PHI Painleve equation. In this approach, the ap- 
pearance of noncompact integrals in the replica limit is a 
mystery. Here, we obtain the resolvent from the solution 
of the replica limit of the Toda lattice equation which 
relates the resolvent to both fermionic and bosonic parti- 
tion functions. The noncompact integrals then arise from 
the bosonic partition function. This method is then gen- 
eralized to the partially quenched case, where the replica 
limit of a generalized Toda lattice equation |29| relates 
the resolvent to a fermionic partition function and a par- 
tition function with both fermions and bosons. 

The partition function of the chUE with n fermionic 
replica flavors with mass m and Nf additional fermionic 
flavors with masses mi, • • • , mjv, is defined by 

N f 

z chUE = ( det n(£, + m )Y[ det(D + m k )). (10) 

k=l 



The Dirac operator is given by 



D 



W 




(11) 



where W is a rectangular I x (I + v) matrix so that D 
has exactly v zero eigenvalues. The average is over the 
probability distribution of the matrix elements of W. We 
will consider the microscopic limit of this partition func- 
tion where the thermodynamic limit is taken for fixed 
x = mUN and x k = m k 'SN. Here, N = 21 + v and 
£ = lim m ^o,JV^oo dm log Z(m)/N. In this case the par- 
tition function reduces to the unitary matrix integral |22] 



Z v n>N M{xi}) = l- [ dUdet'Ue^u+MU*^ (12) 

L-n JU<=U(n+N f ) 



where M = diag(x, • • ■ , x, x\, ■ ■ ■ , xm } ) is the mass matrix 
and an additional normalization factor C n = Y[ k= i{k— 1)! 
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has been included for later convenience. In the micro- 
scopic limit, this partition function coincides with the 
partition function of QCD |30|j2^ ]. The replica limit of 
the resolvent is then given by 

G(x,{ Xi }) = lim -d x logZ» N ,(x,{xi}). (13) 

We first analyze the case of no additional massive flavors. 

The quenched resolvent. In this case the partition func- 
tion is given by 



KM 



1 



Cn 



dUdet»Ue^ u+u ^ 



(14) 



n JUGU(n) 



The unitary matrix integrals are well-known |3^ , ^5| , ^6| , ^7| 



fc=l, 



(15) 



This expression is the r-function of an integrable hierar- 
chy and satisfies the Toda lattice equation [plpS|j3^.|33| 



r A1 2 1 '7V t \ 2 Z n+l( x ) Z n-l( x ) 

[xd x ] \ogZ n (x)=nx [ Z -(xW ' 



(16) 



The choice of the normalization constants C n is consis- 
tent with the large x asymptotic behavior of the l.h.s. of 
( |l6| ) which is given by nx. 

In the replica limit, n — > 0, the partition function can 
be expanded as 



Z v n (x) = l + na" 1 {x)+0{n 2 ). 



(17) 



In this limit, Z^\(x) appears in the r.h.s. of the Toda 
equation, which can be interpreted as the partition func- 
tion with one bosonic replica. With normalization con- 
stant fixed by the large x limit of (|q) it is given by JO] 



Z v _ l {x) = 2K v {x). 



(18) 



For the replica limit of the partition function we thus 
obtain the differential equation 

n [xdx] 2 <j1(x) = 2nx 2 I v (x)K v (x). (19) 

Using the identity 

d x x 2 [K 1/ (x)I v (x)+K v - 1 (x)I v+1 (x)]=2xK 1/ (x)I v (x), (20) 
this differential equation can be integrated resulting in 



al{x) 



dt[t(K v (t)I u (t) + K^xit)!^)) + -] (21) 



in agreement with the known result (|^) for the resolvent! 
The last term is a solution of the homogeneous differ- 
ential equation with integration constants fixed by the 
small x expansion of the partition function. We notice 
that the discontinuity of the resolvent arises through the 
bosonic integral K„-i(x) via the Toda lattice equation. 



In this way, the analyticity in x of the fermionic replica 
expression for the resolvent |8| is avoided. In fyfl , the 
expression for was obtained from the solution of the 
replica limit of the PHI Painleve equation. As a side 
remark, we notice that K v {x) is a solution of the PHI 
equation for n — > — 1 satisfying the boundary conditions 
for n = — 1. The PHI equation thus encodes the analyt- 
icity properties of the bosonic integral K v {x). 

The partially quenched case. The partition function 
( |l2| ) satisfies the recursion relation J29| 

a Zn+i,N f (x, {xk})ZZ_ ltNf (x, {x k }) 
nx ' ; — - r — tt^; (22) 



[Z v n , Nf {x,{x k W 



N f 



x d Xo x k d Xk logZ^ Nf (x , {x k }) 



k=0 

This recursion relation can be derived from a special case 
of the Hirota equation. For degenerate masses x k the sum 
over the derivatives is only over the different masses. The 
normalization of the partition functions is again fixed by 
the large x behavior of the r.h.s. at fixed x k which is 
given by nx. 

Let us consider the replica limit of ( p2| ) for Nf = 1. 
Then the numerator of the l.h.s. is the product of 

2 



■det 



y 



Iv{x) I v {y) 
xl v+1 (x) yl v+ i{y) 



and the supersymmetric partition function ]iT[ 
z -i,i( x ,v) = yK v {x)I v+ i{y) + xK v+ i(x)I v {y) 



(23) 



(24) 



with a bosonic mass equal to x and a fermionic mass 
equal to y. In the replica limit, Z^ 1 (a;, y) reduces to the 
partition function for one flavor and we have 



Z»{x, y) = I v (y)(l + nal{x, y) + 0{n 2 )). 



(25) 



Both sides of the recursion relation are of order n. The 
quantity a\ thus satisfies the linear differential equation 



(xd x + ydy)xd x a\(x, y) = x 
Using ( pcf ) and the identities 



(26) 



(xd x + yd y )[yl l/ (x)l^ +1 (y) - xl v (y)l v+1 (x)] 
= (y 2 - x 2 )I v (x)I u (y), 



I 2 (y)(xd x +yd y 



K u {x) 



VJ J„(v) 

= -[xK v+1 (x)I v (y) + yK v (x)I v+1 {y)]. 



(27) 



this equation can be integrated with solution given by 



t[I u {t)K v {t) + I v+1 {t)K v ^(t) 



v 



KJf) 

I Ay) 



z li(t,y)]dt. 



(28) 
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The resolvent given by d x cr"(x, y) agrees with earlier re- 
sults obtained by the supersymmetric method H]. 

The partially quenched resolvent can be worked out 
analogously for more flavors. For example we have ob- 
tained the resolvent for Nf massless flavors and v = 
from the replica limit of the Toda lattice equation. 

General Solutions. Away from the replica limit we can 
integrate the Toda lattice equation to determine the par- 
tition functions with arbitrary number of fermions and 
bosons but with some degenerate masses. Based on these 
results we conjecture the general form of the partition 
function for arbitrary masses. Up to a normalization con- 
stant we find for pure bosonic partition functions 



det[ 



J- 1 



K, 



v+j-l 



(%i)]i,j=l, 



n i>4 (* 



2 ) 



(29) 



which in the case of degenerate masses reduces to 
Z-n.o(%) — det[Ku+j-i(x)]. This result agrees with the 
integral representation given in [ pl]p5| . The microscopic 
limit of the chUE partition function with n bosons and 
m fermions satisfies the generalized Toda lattice equation 
(^2|) continued to negative n. Iterating this equation we 
conjecture the following result for the partition function 



Z'_ 



,({**}) = 



det[x? Ju+j-i(xi)]i,j=i,..,n+t 



n 



j>i- 



*?)n 



(30) 



The symbol J v {x\) = {— 1 ) v K u (xi) for i = 1, . . . , n and 
J v (xi) = I v {xi) for i = n + 1, . . . , n + m. The structure 
of this expression is in agreement with the general result 
obtained in jj8|. For v = the result for Z v _ n rn was first 
derived by means of the supersymmetric method [j39| . 
The partition function (|3^) is relevant for the calculation 
of higher order spectral correlation functions |^4|,fl0| . 

Conclusions. The two-point function of the Gaus- 
sian Unitary Ensemble and the quenched and partially 
quenched resolvent for the microscopic limit of the chi- 
ral Unitary Ensemble have been obtained by solving the 
replica limit of a (generalized) Toda lattice equation. 
This method provides a direct analytical understanding 
of the appearance of compact and noncompact integrals 
in the expressions for the resolvent. It also explains the 
appearance of a discontinuity in the fermionic replica 
limit of the resolvent. The family of chUE partition func- 
tions with any number of fermions and bosons and with 
arbitrary masses and topological charge also satisfies the 
generalized Toda lattice equation. On this basis we have 
conjectured a form for this partition function that gen- 
eralizes the result for zero topological charge. 
Acknowledgements: We wish to thank P.H. Damgaard 
and M.R. Zirnbauer for useful discussion and suggestions. 
Note added: The expression for Z^_ nm ({xi}) in eq. 
( ^p| ) has now been verified jyj based on a result in |38| . 
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